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Abstract 

The spontaneous flux generation and magnetic field modulation of the critical 

current in a — it Josephson junction are calculated for different ratios of the 

junction length L to the Josephsen penetration depth Aj, and different ratios 

of the 0-junction length to the 7r-junction length. These calculations apply to 

a Pb-YBCO c-axis oriented junction with one YBCO twin boundary, as well 

as other experimental systems. Measurements of such a junction can provide 

information on the nature of the c-axis Josephson coupling and the symmetry 

of the order parameter in YBCO. We find spontaneous flux even for very 

short symmetric — it junctions, but asymmetric junctions have qualitatively 

different behavior. 
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I. INTRODUCTION 



Although many measurements now support the idea that the gap in the high-temperature 
cuprate superconductors has dominant d x i_ y i character ffl], an orthorhombic superconductor 
such as YBCI2CU3O7-S (YBCO) cannot be described as a purely tetragonal cf-wave. In the 
orthorhombic symmetry group the d x 2_ y2 — and s— wave basis function belong to the same 
irreducible representation, so that one expects that the gap can be described as d x 2_ y 2 — wave 
with some s— wave admixture, d + as. The deviation from purely ci-wave behavior is sup- 
ported by the observations of Josephson pair tunneling from Pb into c-axis oriented 
YBCO, although the variation of the Josephson tunneling strength with twin density, the 
relative areas of the twins, and the possible role of Josephson screening currents, remain 
open questions [|J. It has also been proposed that these results can be explained by the 
development of a complex order parameter at twin boundaries ||, or mixing of states with 
s and d x 2_ y 2 symmetry at the Pb-YBCO interface ||. Part of the problem of interpreting 
these experiments is the complex geometry presented by highly twinned materials. Here 
we consider the simpler case of 8b C~£IX1S YBCO-Pb tunnel junction in which the YBCO has 
only one twin boundary in the junction. Experimental efforts on this system are currently 
underway |7|]. 

In the idealized geometry we describe below (Fig. 1), the Pb/YBCO junction which 
contains a single twin boundary can be described as a "0 — 7r" junction, where the pair 
transfer integral has a relative pair phase of in the left part of the junction (length Lq) and 
it in the right part of the junction (length L n ). Here we are interested in understanding the 
behavior of such "0 — n" Josephson junctions. While we describe our geometry in terms of 
c-axis tunneling in YBCO-Pb junctions, our results apply equally well to the — it a-b plane 
YBCO-Pb planar junctions || and grain boundary junctions [|9|jll| used for phase sensitive 
tests of the symmetry of the high-T c superconductors. In particular, while traditionally one 
considers the magnetic interference pattern, we will also be interested in spontaneous flux, 
flux generated by the self-screening Josephson currents in the absence of a drive current 
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and an external field H e , which can be directly probed by, for example, a scanning SQUID 



microscope ||. As noted in Ref. ||TT|, the two measurements are complementary, since while 
the magnetic field interference patterns are only strikingly different between a — it junction 
and a conventional — junction in the short junction limit L ~ Aj, appreciable spontaneous 
flux only appears in the — tx junction in the long junction limit L ~ Xj. 

The maximum Josephson current J c which a junction can carry versus an external mag- 
netic field H e applied parallel to the plane of the junction is called the "magnetic interference 
pattern". The interference pattern for a short (L <C Aj) — it junction has been discussed 
by Wollman et al. ||. Xu et al. [|TTJ have estimated the crossover to the limiting case of a 
long (L ^> Aj) — 7r junction. With the approximations made by Xu et al., the magnetically 
modulated critical current of a long — ir junction is basically identical to the conventional 
— junction due to entrance of a half-flux quantum vortex, and the solution which con- 
tains spontaneously generated flux ceases to be a global minimum of the free energy for 
short junctions. Our exact numerical solutions, however, show that there is still a "dip" 
in the center of the diffraction patterns even for junctions as long as 10Aj, and that the 
symmetric junction should contain spontaneously generated magnetic flux even for L < Xj. 
Here we are primarily interested in determining the behavior of a — 7r junction for inter- 
mediate values of L/Xj and different L Q /L n asymmetries. We believe that the dependence 
on L/Xj of both the field- modulated critical current and the spontaneous flux generation 
of a Pb-YBCO junction with one twin boundary can provide important information on the 
c-axis coupling. 

In the following, we discuss a numerical method for calculating the spontaneous flux 
generation and magnetic interference patterns for arbitrary junction lengths L/Xj and junc- 
tion asymmetries L Q /L n . We then compare results for the critical current I c versus external 
magnetic field H e for the case of a traditional — junction with the symmetric — it case as 
a function of the reduced junction length L/Xj. Following this, we focus on the Josephson 
screening currents by examining the flux generated by these self-currents in the absence of 
an external current and magnetic field. As previously noted, this self generated flux can 



be probed using a scanning SQUID microscope |J, providing important information on the 
nature of the c-axis Josephson coupling. We then turn to the interesting case in which the 
areas of the two twin regions differ and again look at the self-flux versus L/Xj. There is a 
qualitative difference between the results for the spontaneous flux generation which occurs in 
asymmetric junctions compared to symmetric junctions. We show briefly how these results 
can be extended to a junction with multiple twin boundaries. We conclude by discussing 
what this type of experiment can tell us about the Pb-YBCO c-axis Josephson tunneling. 



II. THE CALCULATION 



The geometry and gap structure which we envision is illustrated in Fig. 1. In this 



idealized geometry the upper YBCO strip contains a twin boundary which separates 
the YBCO-Pb junction into two regions. As schematically illustrated in Fig. 1, we have 
taken the phase of the larger lobe of the d + as gap in region < x < Lq to be positive. 
Then, if the phase is locked across the grain boundary as one expects [[Uj, the larger lobe in 
region L < x < L will be negative. Thus the pair transfer integral between the YBCO and 
the Pb leads to a relative Josephson pair phase of for < x < L and tt for L < x < L. 
Following the notation of Owen and Scalapino JL2], we describe a junction (Fig. 1) with 



width w small compared to Aj in the y-direction, of length L in the x-direction, carrying 
a total current Itot in the x-direction, in an external magnetic field H e oriented in the y- 
direction. We include the effects of the and tt junctions by introducing an extra phase 
angle 6(x) which takes the value or tt. The superconducting phase difference across the 
junction is then just the solution of the Sine-Gordon equation: 

d 2 s i 



dx 2 Xj 



sm((f)(x) + e(x)). (2.1) 



Applying Ampere's Law with a contour of integration around the perimeter of the junction 
in the xy plane leads to the boundary condition 

H(L) - H(0) = 4nI tot /cw. (2.2) 



Contours of integration in the yz plane circling the leads at x = and x = L lead to 



H(L) + H(0) = 2H e . (2.3) 



The Josephson penetration depth is given by 



(2 - 4) 



where d is the sum of the Pb and YBCO \ a b penetration depths plus the thickness of the 
insulator layer between the two superconductors, and j c is the Josephson critical current 
density The current per unit length through the junction is given by: 

j(x) = wj c sin(<f>(x) +9(x)). (2.5) 

The gradient of the phase is in turn related to the field H in the junction by: 

„. s he 86 

H(x) = — ~. 2.6 
If we redefine the parameters H(x) and I to t as 

h(x) = —H(x) = Xj-, (2.7) 

and 

Hot = (2.8) 

the boundary conditions become 

h(L) = i tot + h e (2.9) 

h(0) = -i m + h e . (2.10) 



Defining a dimensionless coupling parameter 



the differential equation turns into a difference equation on a grid of size Ax: 

<p n+1 - 2<p n + n _! = - sin (0 n + n ). (2.12) 

a 

The boundary conditions are then described as difference equations, where rij is the total 
number of junctions: 

<t> nj ~ <t> nj -i = Hot ~^ ke (2.13) 



itot + h e ^ 14 ^ 



'a 

These coupled difference equations are solved using a relaxation method to find the solution 
4>(x). The free energy of this solution is given by: 

<«*> + '<*»4(i)> (2 - i5 » 

Written as a difference equation, the free energy for the vortex solution becomes: 



„ hj c wAx n 
F v = 



2e E (i - COS (<^ + °n) + |(0n+l " 0n) 2 ) (2-16) 

while the free energy of the no- vortex (0 = everywhere) solution is: 

Fo = M^"£ (1 _ cosf y (2.17) 

To obtain the critical current, the current I to t through the junction is increased in steps, 
iterating the solution until either the convergence criterion 

6 > {[(0i -<h-(i- K)/^ 2 ) 2 + (0„ 3 - n ._ 1 - (i + h^/a 1 ' 2 ) 2 

rij—l 

+ E (^n+i + 0n-i - 20„ - s^n(0 n + 6 n ) / a) 2 } / n^ 2 / IT (2.18) 

n=2 

is satisfied, with e normally chosen to be 1CT 6 , or until the solution diverges with | (f>(rij/2) \> 
25. For particular values of the parameters, we compare solutions with e taken to be 10~ 6 , 
10~ 7 , and 10 -8 to check the convergence of the solution. The largest value of Itot for which 
the iterations converge is taken to be the critical current J c . For the numerical calculations 
shown in the next section, we took L = 10, Ax = 0.1, rij = Lj Ax = 100, which meant 
that the dimensionless coupling constant ranged from a = 10 4 for L/Xj = 1 to a = 10 2 for 
L/\j = 10: Ax was always much smaller than L , L n , or Aj. 
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III. RESULTS 



Fig. 2 compares our results for the magnetic interference pattern for a — junction 
(a) and a symmetric — tt junction (b) for various values of the reduced length L/Xj. Here 
we have plotted I c /h versus $/$o? where Ji = j c wL, $ = H e dL, and $o — hc/2e is the 
superconducting flux quantum. In the short junction limit L <C Xj, the critical current can 
be written as: 

IJB) = -± max / sin * + 0(g) + V)<*c , 3.1 

where 6(x) can be or tt, and the maximum value of the expression is obtained by varying 
< ip < 2tt for each value of H e . For the — junction (0(x) = 0) this expression reduces 
to: 

Ic = I 3zn(7r$/$ ) | 

while for the symmetric — it junction (9(x) = 0, x < L/2; 6(x) = tt,x > L/2) this becomes: 

I c _ g zn 2 (7r$/2$ ) 

Ji ~ | 7T$/2$ | • ^ 
These equations are plotted as the solid lines in Figure 2. Also included in Figure 2a is the 
result of Owen and Scalapino for the reduced length L/Xj = 10. For the — junction the 
effect of increasing L is to reduce the height of the central peak in the magnetic interference 
pattern, and to reduce the amplitude of the successive oscillations as H e is increased. For 
the — tt junction, increasing L tends to reduce the depth of the minimum at H e = 0, and 
as L/Xj — > oo the curves for the — junction and the — tt junction become identical 



as discussed by Xu et al. |TTJ. However, even for L/Xj = 10, we find, as shown in Fig. 



2(b), that the critical current initially increases with flux until $/$o — 0.5. This differs 



from the L/Xj = 10 result shown in Ref. ||11||. We believe that this difference results from 



the approximation made by Xu et al. that the presence of a tt vortex simply changes the 
phase in the it section of the junction by tt, effectively changing it into a 0-junction. This 
approximation is valid on length scales large compared to Xj, but, as shown in Fig. 2b, leads 



to an experimentally detectable difference in the solutions even for L as long as 10A j. For 
short and medium length junctions, Xu et al. took the analytical long-junction expression 
and made the approximation that the effect of shortening the junction would be simply to cut 
off this expression. However, this is not a valid solution of the modified sine-Gordon equation 
and the boundary conditions. It also leads, as we show below, to dramatic differences in 
both the free energy and the total flux for the 7r- vortex solutions for short- to medium-length 
junctions. 

Fig. 3 addresses the question of spontaneous flux generation in — n junctions as a 
function of reduced length L/\j and the degree of asymmetry (L n /L). Fig. 3a plots the 
ratio of the free energy of the state with some spontaneous flux to the state with no flux. 
The dashed line is the result of Xu et al, using the approximation outlined above. The 
solid lines are our full results. Note that, in contrast to the results of Xu et al, for the 
symmetric — it junction (L n /L = 0.5) the state with spontaneous flux always has lower 
energy than the state with no flux, and thus some self generated flux should therefore be 
present for all values of L/\j. Fig 3b plots the spontaneous flux vs. L n /\j. The inset in 
this figure shows that the spontaneous flux for a symmetric junction follows a power law 
dependence on junction length for short junctions. In this limit the phase has only small 
deviations from an average value of n/2, and the spontaneously generated field increases 
linearly towards the center of the junction from a value of zero at the edges. A simple 
geometrical argument, expanding the Sine-Gordon equation about = tt/2, implies that for 
short junctions the spontaneous flux should be given by &/&o — (L/ Aj) 2 /8ir. This relation 
is plotted as a straight line in the inset of Fig. 3b. However, for the asymmetric case, the 
state with no flux has the lowest free energy for short junctions, and as shown in Fig. 3b, 
there is no spontaneously generated flux, up to a critical value of L n /Xj ~ 1. The amount 
of spontaneously generated flux approaches $o/2 as L n gets larger, and rate of increase of 
flux at the onset of flux generation increases for less symmetrical junctions. 

The onset of spontaneous flux generation is also apparent in the magnetic interference 
patterns, as shown in Fig. 4. Here we hold the asymmetry factor L n / L fixed at 0.25, and vary 

8 



the length of the junction L n / Xj. For small L n / Xj there is no spontaneous flux generated, 
and the magnetic interference has a minimum at zero field. However, as L n / Xj approaches 
unity, there is an abrupt shift to a magnetic interference pattern with a maximum at zero 
field. The solid line in Fig. 4 is the short-j unction limit calculated from Eq. 3.1, which 
is a good indicator of the actual interference pattern only until flux generation becomes 
energetically favorable. 

The behavior of junctions with multiple twin boundaries can be calculated using the 
same techniques, as long as the junction width w is small compared to Xj. We show the 
results of one such set of calculations in Fig. 5. Here we have assumed that there are 
100 twins in the junction, with the probability of a particular twin having 0— or 7r— phase 
randomly distributed with the probability of 7r being 0.25. The circles in Fig. 5a show the 
randomly distributed set of 0— and tt phases used in this calculation. The solid lines in Fig. 
5a show the solution for the phase angle <fi at H e = and I = I c , for various values of L/Xj. 
These curves show, as expected, that the gradient of 0, and therefore the flux penetration 
into the junction, is spread throughout the junction for small L/Xj, but is localised at the 
junction edges for L/Xj ^> 1. The solid line in Fig. 5b shows the predicted interference 
pattern from Eq. 3.1 in the short-j unction limit. The other symbols show the result of the 
full numerical solution for L/Xj = 2,4, and 10. Even for this small number of randomly 
distributed twins, the interference pattern, although reduced in critical current, is similar to 
that for a single 0-junction (Fig. 2a) especially for larger values of L/Xj. This figure shows 
how important it is that complete knowledge of the distribution of twins be available for the 
correct interpretation of experimental work on c-axis tunneling from twinned samples. 

IV. CONCLUSION 

As discussed by Sun et al. ||, the results of their Pb-YBCO c-axis Josephson tunneling 
experiments can be interpreted as showing that the order parameter in YBCO must have 
some s-wave component. While this could be consistent with an order parameter having 
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d x 2_ y 2 + as symmetry in an orthorhombic material, the interpretation of the experiments 
done on twinned materials depend on details of the multi-twinned geometry. Here we have 
analyzed a simpler situation involving a single twin boundary. In this case, if the d x i_ y i 
contribution is dominant, giving the type of order parameters illustrated in Fig. 1, one 
should observe the behavior we have discussed. In particular, if the junction areas for the 
two twin regions are equal, one has the conventional minimum in the critical current at 
zero magnetic field as shown in Fig. 2(b), if L/\j is sufficiently small, and provided the 
current density is uniform on a scale set by Xj. For the asymmetric case in which the areas 
of the two twin regions differ, more striking behavior should appear. There should be a 
sudden shift of the interference pattern from one with a minimum at zero field, to one with 
a maximum (Fig. 4), as L n /Xj goes through unity. One way to vary L/\j would be by 
changing Aj with temperature. Another would be to apply a magnetic field H x to vary 
\j(H x ) = [hc 2 /(^edI(H x )/wL)] l l\ with I{H X ) = h \ sin{>K<$> x /<$> ) \ / \ ir$ x /$ | , and 
§ x = H x dw. Alternatively, or in addition, one could use a scanning SQUID microscope to 
measure the spontaneous flux generated by the self-screening Josephson currents. For the 
case of equal areas, this spontaneous flux rises continuously as L w / A j increases. However, for 
nonequal areas one finds no spontaneous flux generation below a critical value of L n / Aj ~ 1 
and then a sharp increase in flux should occur. We believe that if these types of behavior 
are observed, they would show that the Pb-YBCO c-axis Josephson coupling is consistent 
with a description of the high T c superconductors as having dominant d x 2_ y 2 symmetry. 
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FIGURES 

FIG. 1. (a) Junction geometry showing the directions of the external current flow I tot and mag- 
netic field H e . The dashed line marks the twin boundary, and the gap functions are schematically 
illustrated such that the basic Josephson pair phase arising from the pair wave function overlap is 
for < x < Lq and it for Lq < x < L. (b) Relative pair tunneling phase across the junction as a 
function of position. 

FIG. 2. The dependence of the junction critical current I c /h as a function of the field in the 
junction $/$o> for a series of junction lengths L/Xj, for a — junction (a) and a symmetric 
— 7T junction (b). The solid curve is the analytical result in the short junction L <C Xj limit. 
The dashed curve in (a) is the result for L/Xj = 10 from Ref. [12]. Curves for successive values of 
L/Xj have been offset vertically by 0.5 units for clarity. 

FIG. 3. The ratio of the free energy for the solution with spontaneous flux F v to that with 
no flux Fq (a), and the spontaneous total flux in the junction $/$o (b), for a — 7r junction as 
a function of the length of the tt junction L w /X,j, for different asymmetries L w /L. The inset of 
(b) shows the results for a symmetric junction on a log-log scale. The straight line is the relation 
$/$o = (L/Xj) 2 /8tt. 

FIG. 4. The dependence of the critical current I c /h as a function of field in the junction 3>/3>o> 
for a junction with asymmetry L w /L = 0.25, for three different junction lengths: L T /X j = 0.5, well 
below the critical length for spontaneous flux generation; L w /Xj = 0.75, near the critical length; 
and L n /Xj = 1, above the critical length. The solid curve is calculated using the short junction 
approximation, Eq. (3.1). Curves for successive values of L/Xj have been offset vertically by 0.5 
units for clarity. 
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FIG. 5. (a) The circles show 100 randomly distributed and it phases, assuming a probability 
of 7r-phase=0.25. The solid lines show the solutions for cf) at I = I c and H e = 0. (b) The solid 
curve is the magnetic interference pattern for a junction with this distribution of phases in the 
short junction limit, Eq. (3.1). The solid points are the result of the full numerical calculation for 
various values of L/Xj. Successive L/Xj curves have been offset vertically by 0.2 units for clarity. 



13 



